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Introduction 


V 

In  1939  Flory  announced  the  principle  of  equal  reactivity  of  polymer 
chains  (1).  The  principle  maintains  that  the  ability  of  a  chain  molecule 
to  form  larger  molecules  is  independent  of  the  chain  length.  The  argtmient 
was  very  significant  at  that  time  for  two  reasons:  First,  the  existence  of 
high  molecular  weight  polymers  was  considered  by  some  to  be  improbable. 

These  proponents  argued  that  the  high  molecular  weight  species  must  have  a 
higher  diffusion  constant,  relative  to  the  lower  molecular  weight  species, 
which  would  eventually  limit  the  extent  of  polymerization.  Studies  on 
condensation- type  polymers  with  random  coils  demonstrated  the  general 
validity  of  the  principle  of  equal  reactivity  (2).  The  second  signifi¬ 
cant  result  was  that  the  distribution  of  molecular  weights  could  be 
readily  calculated  from  kinetic  parameters  that  were  independent  of  chain 
length. 

The  major  significant  deviations  from  this  principle  seem  to  occur  only 
with  very  low  molecular  weight  species  where  inductive  or  polarization  effects 
can  be  transmitted  across  a  small  molecule.  Indeed,  Lenz  (3)  cites  a  number 
of  cases  where  the  reactivity  of  a  functional  group  on  a  difunctional 
monomer  depends  on  whether  or  not  the  other  functional  group  has  reacted. 
Challa  (4)  showed  a  similar  effect  on  the  rate  of  polymerization  of  poly- 
(ethylene  terephthalate)  where  the  reactivity  of  the  "monomer"  bis- 
(B-hydroxyethyl)  terephthalate  differs  from  that  of  higher  oligomers  and 
polymer.  Peebles  and  Wagner  (5)  examined  the  transesterification  kinetics 
of  dimethyl  terephtalate  and  ethylene  glycol  in  the  presence  of  a  catalyst. 
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They  concluded  that  the  second  transesterlflcatlon  reaction  was  some  three 
times  faster  than  the  first  transesterification  reaction’.  The  overall 
transesterification  of  the  terephthalate  monomer  must  somehow  involve  some 
type  of  molecular  transition  state  containing  the  catalyst  because  attempts 
to  verify  the  faster  transesterification  reaction  with  the  partly  transester- 
ified  material  resulted  in  absolutely  no  reaction  (6).  Thus,  the  documented 
"violations"  of  the  principle  of  equal  reactivity  seems  to  be  restricted  to 
very  small  chain  segments  of  random  coil  molecules. 


Several  calculations  have  been  made  to  determine  how  departures  from 
the  principle  of  equal  reactivity  would  influence  the  molecular  weight  dis¬ 
tribution  of  a  polymer.  Nanda  and  Jain  (7)  assumed  that  the  second  order 
reaction  constant  for  a  condensation  polymer  was  proportional  to  the  molecular 
weight  of  the  formed  polymer  through  the  relation  k(l+  where  k  and  \ 
are  constants  and  i^  is  the  number  of  monomers  in  a  chain  of  length  i,, 

to  obtain  an  analytical  solution  for  the  molecular  weight  dis¬ 
tribution.  Magat  (8)  developed  simultaneous  integral  equations  for  stepwise 
polymerizations  without  termination  for  a  number  of  conditions  where  the  rate 
constant  for  propagation  varied  with  chain  length.  His  approaches  are 
limited,  however,  in  that  the  precise  polymerization  conditions  must  be 
specified  before  the  equations  can  be  evaluated. 


A  new  series  of  stiff  rod-like  polymers  (9)  have  recently  been  synthesized 
and  characterized  with  a  repeating  unit  of  the  form 


ICE' 


or 


where  the  "X"  atom  can  be  oxygen,  sulfur,  or  singly  hydrated  nitrogen. 


The  condensation  of  the  functional  groups 


appearsto  occur  simultaneously  to  form  the  benzobisoxazole  because  intermediate 
functionalities  have  not  been  detected  by  chemical  or  spectral  means  (9e}. 

These  observations  suggest  that  the  interacting  molecules  must  be  co1 linear 
(or  nearly  so)  before  reaction  can  take  place.  Furthermore,  once  reaction 
has  commenced,  the  entire  series  of  reactions  to  form  the  cyclic  structure 
must  occur  quite  rapidly.  If  such  is  the  case,  then  the  rate  determining 
step  is  the  ability  of  rods  to  rotate  into  a  col  linear  configuration  rather 
than  the  coupling  step.  Here  the  principle  of  equal  reactivity  is  not  violated; 
the  reactivity  of  each  end  is  independent  of  the  length  of  the  chain  to  which 
it  is  attached.  The  purpose  of  this  paper  is  to  calculate  the  molecular  weight 
distribution  for  randomly  oriented  rod-like  molecules  which  must  diffuse  or 
rotate  into  a  col linear  configuration  prior  to  the  coupling  reaction. 
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Mathematical  Development 


The  reaction  between  molecules  of  size  j[  reacting  with  molecules  of  size 
i  to  form  a  molecule  of  size  i  +  i  irreversibly  and  without  by-product 

(1) 

is  postulated  to  have  a  rate  of  formation  proportional  to  k  where  k  is 

the  specific  second  order  rate  constant  for  coupling,  Cy^  is  a  function  of  the 
molecular  lengths  j.  and  j.  and  is  a  measure  of  the  ability  of  molecules  to  become 
coll  inear,  and  the  concentrations  of  molecules  of  size  i_  and  i- 


To  simplify  the  notation,  the  constant  k  is  adsorbed  into  the  time  derivative. 
Thus  from  eq  (1),  the  rate  of  formation  of  x^  is  given  by 


'j.i-jYi-j  -  "I  jf,  'u  ''j 


The  monomer  x^  is  considered  to  react  at  the  rate 

CB 

*1  ■  -  ’'i/.,  '14  "j  (3) 

These  rate  formulas  can  be  summed  to  give  the  following  time  derivative  functions 


-  0 


since  no  by-product  is  formed. 


00  2  •  • 

^  i  x^  ■  s  i  x^  r  J  c,,  X, 
i-1  ^  i-1  ^  J-1  'J  J 


If  the  variable  c^can  be  specified,  then  the  molecular  weight  distribution 
and  the  moments  of  the  distribution  can  In  principal  be  determined  by 
computer  calculations. 

Because  many  polymers  have* degrees  of  polmerizatlon  In  excess  of  100, 
and  perhaps  Into  the  thousands  or  millions.  It  Is  not  reasonable  to  solve 
the  excessive  number  of  simultaneous  differential  equations  given  by  (2) 
and  (3).  However,  If  some  type  of  scaling  law  can  be  found  which  will  permit 
calculation  of  the  distribution  for  a  low  number  average  molecular  weight 
distributions,  then  computer  calculations  are  feasible. 

The  generalized  exponential  distribution  function 

F(1,m.k)  =  my''^"’  f''’‘'Cexp(-yl'")]/p(k/m)  (7)  Ci 

W(1,m,k)  •  l’^[exp(-yl'")]/p[(k+l )/m]  (8) 

where  <i>^  *r  [(k+1  )/m]/y^^'"p(k/m), 

can  be  rewritten  In  a  scaling  law  form  >w1th  as  the  Independent  variable. 

[  i  ^  «p  [.  (,1/<1>„r]  (  9) 

with  g  ■  {pC(k+l)/m]}/f(k/m) 

This  distribution  Is  a  generalized  form  of  the  Schulz  distribution  (m  ■  l)f 
the  Tung-Welbull  distribution  (k  ■  m-1,  k  >  0}»  and  a  good  approximation  of 
the  most  probable  distribution  Cm  •  k  •  1)  at  large  values  of  <1>j^. 

The  next  question  to  consider  Is  the  form  of  the  function  c^j  which 
determines  the  ability  of  the  molecules  to  rotate  Into  a  colllnear  conformation. 
Doe  and  Edwards  (9)  suggest  that  the  rotational  diffusion  constant  of  a  rod 


V 
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depends  on  (In  1)/i  in  dilute  solution  and  higher  Inverse  powers  of  1  In  more 

concentrated  solutions,  whereas  Bird  (10)  suggests  that  the  rotation  of  a 

rigid  dund)be11  depends  on  the  square  of  the  length  (1  ).  The  assumption  can 

then  be  made  that  If  rod  1^  were  held  fixed  In  space  while  rod  j.was  allowed  to 

•N 

rotate  freely,  then  the  rate  of  rotation  would  be  proportional  to  j  , 

where  N  Is  an  adjustable  parameter.  If  rod  iwere  held  fixed  and  rod  1,  allowed 

•N 

to  rotate,  the  rate  Is  proportional  to  1"  .  The  overall  average  rate  Is  then 
c^j  *  (r*^  +  j‘”)/2 

This  equation  can  then  be  substituted  Into  the  set  of  simultaneous  differential 

equations  2r^-^.  The  initial  concentration  of  monomer,  was  selected  to  be 

unity.  Initially,  was  taken  as  the  Independent  variable.  A  set  of  K 

simultaneous  differential  equations  ’^tnd 

concentration  with  time  of  the  Independent  variable,  was  set  up  for  solution 

by  a  computer.  The  value  of  K  was  selected  to  be  large  enough  so  that^1x^).999 

over  the  range  of  Integration  considered.  A  Runga-Kutta  method  was  used  to 

obtain  values  of  xj,  from  1_  *  Ind  to  1_ *  K .  After  each  Integration,  an  Interger 

value  of  L  =  0.8K  was  tested  to  see  whether  LXl<0.0005.  If  the  test  failed, 

a  larger  value  of  K  was  selected,  the  program  retreated  to  the  results  obtained 

on  a  previously  performed  Integration,  and  computation  continued.  This 

procedure  allowed  only  a  small  ninnber  of  equations  to  be  solved  at  the  beginning 

of  the  solution,  but  as  <i>^  increased,  more  equations  were  solved,  and  hence 

more  computer  time  required.  Also,  as  <1>j^  increased,  values  of  ^nd  decreased 

•  « 

and  eventually  became  too  small  to  permit  reliable  values  of  to  be  com¬ 
puted  as  evidenced  by  fluctuations  in  As  x^^^  approached  this  limiting 

value,  a  new  Independent  variable  was  selected  and  concentrations  of  all  1, 
variables  below  were  arbitrarily  set  to  zero.  Extensive  testing  of  this 

procedure  by  choosing  various  values  of  1_  for  the  Independent  variable  exhibited 
no  effect  on  the  calculated  distribution  provided  that 
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Results 


Integration  of  eq's  2,  3,  5  and  6  were  performed  for  integer  values 
of  N  in  eq  li  as  a  function  of  conversion.  Values  of  <i>j^  varied  from  unity 
(initial  condition)  to  <i>jj  values  in  excess  of  50.  From  the  X.  values  so 
computed,  values  of  x^<i>^/rxj  were  obtained.  For  all  values  of  N  used, 
except  for  N  «  0  (the  most  probable  distribution)  the  scaling  law  held  to 
within  the  precision  of  the  data  when  <i>|^  values  exceeded  ten.  The 
results  are  presented  in  Fig.'s  1  and  2  for  the  distribution  as  a  function 
of  i/^i^n  aod  in  Fig.  3  for  the  heterogenity  index  *  function 

of  conversion. 


Discussion 

The  empirical  observation  that  the  scaling  laws  F(1}  <i>^  and  H(1}<1>^ 

are  functions  only  of  V<1>„  means  that  Inexpensive  computer  calculations 

can  be  performed  to  evaluate  these  functions  at  small  values  of  <1>.. 

then  " 

The  function  can^be  used  directly  to  evaluate  the  molecular  weight  distribu¬ 
tions  for  large  values  of  <1>„.  This  technique  can  be  used  for  other  schemes 
relating  kinetics  to  reactor  designs*  etc.  where  numerical  calculations  are 
required. 

When  the  value  of  N  In  equation  11  Is  nonzero,  the  rate  of  polymerization 

eq  5,  decreases  with  Increasing  conversion.  As  N  Increases,  the  rate  of 
polymerization  decreases  at  a  faster  rate.  Both  the  scaled  mole  fraction  and 
weight  fraction  distributions  show  a  maxlmun  in  the  range  1/<1>n  of  0.6-0. 7, 
the  maximum  shifting  to  higher  values  with  Increasing  N.  The  breadth  of  the 
distribution,  however,  becomes  narrower,  and  reaches  an  asymptotic  value  at 
<1>j,  values  of  approximately  10.  The  application  of  these  calculations  to 
the  observed  distributions  of  Poly([benzo(1,2-d:5,4-d')  blsoxazole-2,6  -dlyl] 

-1, 4-phenyl ene)  has  been  discussed  by  Cotts  and  Berry  (9e). 
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Legend  for  Figures 

Figure  3.  The  scaled  "lole  fraction  of  polymer  as  a  function  of  for 

various  values  of  N  in  eq  n 

Figure  2.  The  scaled  weight  fraction  of  polymer  as  a  function  of  V<i>^  for 
various  values  of  N  in  eq  11 

Figure  3.  The  heterogeneity  index  ®  function  of  <i>j^  and  N. 
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